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CRITICAL  PHENOMENA  ON  SURFACES: 
RECONSTRUCTIONS  AND  CHEMISORBED  LAYERS 

W.N.  Unertl 


ABSTRACT 

\ 


This  article  reviews  recent  experimental  studies  of  the  phase  transitions  which 
occur  on  reconstructed  surfaces  and  in  layers  of  chemisorbed  atoms.  Emphasis  is 
given  to  systems  which  exhibit  continuous  transitions.  The  experimental  results  are 
presented  in  the  context  of  current  theories  which  emphasize  the  roles  of  symmetry, 
dimensionality  and  universality.  The  limitations  currently  imposed  by  sample  quality 
and  experimental  techniques  are  also  discussed.  Results  from  individual  systems  are 
surveyed  and  promising  future  directions  including  studies  of  finite  size  effects  are 
discussed. 
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1.  INTRODUCTION 

The  concepts  of  symmetry,  dimensionality  and  universality  are  central  to 
contemporary  theories  of  critical  behavior.  Indeed,  some  classes  of  transitions  are 
predicted  to  occur  only  in  two  dimensional  systems.  This  fact  has  led  to  an  outburst 
of  experimental  activity  to  study  critical  phenomena  on  real  surfaces.  Second-order 
phase  transitions  have  special  importance  in  this  context  and  the  purpose  of  this 
review  is  to  provide  a  survey  of  the  extent  to  which  recent  experiments  support 
current  theoretical  descriptions  of  second-order  phase  transitions. 

Surface  phase  transitions  are  also  of  interest  because  they  can  effect  such 
diverse  phenomena  as  heterogeneous  catalysis  and  thin  film  growth.  Surface  diffusion 
[1],  adsorption  and  desorption  rates  [23,  chemical  reactions  on  surfaces  [3],  and  the 
degree  and  type  of  order  in  thin  films  [43  are  other  examples  of  processes  which  can  be 
Influenced  by  phase  transitions. 

The  equilibrium  properties  of  a  surface  system  are  conveniently  summarized  on  a 
phase  diagram.  This  article  is  concerned  only  with  closed  systems;  i.  e.  those  for 
which  the  number  of  particles  in  the  surface  layer  remains  constant.  The  phases  for 
closed  systems  are  usually  represented  on  a  Temperature  vs.  Coverage  (T  vs  0)  phase 
diagram.  Figure  1  is  a  typical  example,  in  this  case  for  selenium  adsorbed  on  the 
Ni(  1 00)  surface  [51.  The  ranges  of  T  and  0  over  which  the  various  phases  exist  are 
delineated  by  phase  boundaries  which  can  be  either  first-  or  second-order.  The 
crystallographic  structure  of  each  phase  and  the  precise  location  of  its  phase 
boundaries  are  determined  by  the  interactions  of  adsorbate  atoms  with  each  other  and 
with  substrate  atoms.  Figures  2a  and2b  show  the  atomic  structure  of  the  simple  p(2  x 
2)  and  c(2  x  2)  phases  of  Figure  I.  The  sequence  in  which  these  phases  are  formed  as  0 
increases  provides  direct  information  about  the  relative  interaction  strengths 
between  Se  atoms  as  shown  in  Figure  2c.  For  example,  in  the  perfect  p(2  x  2)  phase, 
nearest-neighbor  and  next-nearest-neighbor  adsorption  sites  are  not  occupied  by  5e 
atoms.  Thus,  the  interaction  E3  must  be  more  attractive  than  either  Ej  or  E2.  Since 
the  c(2  x  2)  phase  forms  next,  E2  must  be  more  attractive  than  Ej.  Therefore,  we 
conclude  that  E3  <  E2  <  Ej.  The  numerical  values  of  these  energies  can  be  estimated  by 
fitting  calculated  phase  boundaries  to  the  experimental  boundaries  using  various 
statistical  mechanical  techniques  [53.  Such  fits  provide  one  of  the  few  methods 
currently  available  for  determining  the  interaction  energies  between  adsorbed  atoms. 
This  is  one  of  the  major  motivations  for  studying  surface  phase  diagrams. 

The  variety  of  possible  surface  phase  transitions  is  immense.  This  article  is 
restricted  to  only  a  small  subset  of  these:  ( I )  phase  transitions  which  occur  on  clean, 
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reconstructed  surfaces  (e.g.  Au(  It  0)  (1x2)  and  W(  100)  c(2x2)),  and  (2)  phase 
transitions  in  ordered  layers  of  chemically  adsorbed  atoms.  Other  important  systems 
include  physisorbed  layers  [6],  intercalation  layers  [7],  liquid  crystals  [8]  and  lipid 
films  [91  Furthermore,  we  discuss  only  systems  for  which  the  phase  diagram  and  its 
phase  boundaries  have  been  experimentally  determined.  Finally,  only  adsorption 
systems  which  can  reasonably  be  described  by  lattice  gas  models  are  included;  1.  e. 
those  for  which  adsorption  is  at  specific  sites  on  the  surface.  Thus,  commensurate- 
incommensurate  transitions  and  sequences  of  compressed  phases  are  outside  the  scope 
of  this  article  [10]. 

The  experiments  discussed  here  rely  entirely  on  Low  Energy  Electron  Diffraction 
(LEED)  to  characterize  the  degree  of  order.  The  reader  is  assumed  to  be  familiar  with 
the  basic  ideas  of  diffraction  theory  and  the  reciprocal  lattice  [1 1,12]  and  to  have  a 
working  knowledge  of  the  notation  used  to  describe  two  dimensional,  periodic 
structures  [13]. 

2.  THEORETICAL  BACKGROUND 

This  section  summarizes  some  of  the  basic  theoretical  background  needed  to 
place  the  experimental  studies  of  continuous  order-disorder  phase  transitions  into  the 
context  of  modern  theories  of  critical  behavior.  Since  this  presentation  must  be  brief, 
references  are  given  throughout  to  more  detailed  treatments. 

2. 1  Order  Parameter.  Pair  Correlations  and  Correlation  Length 

Three  measurable  quantities  whose  behaviors  in  the  vicinity  of  a  phase  transition 
are  predicted  by  theory  are  the  order  parameter,  the  pair  correlation  function  and  the 
correlation  length.  Consider,  for  example,  the  case  of  a  lattice  gas  of  particles 
adsorbed  on  a  lattice  with  unit  mesh  vectors  a  j  and  a 2  and  lattice  sites  at  rj  -  m  |  ,a  j 
♦  m2ja2  where  m  j  j  and  m2j  are  integers.  At  any  instant  the  probability  nj  that  a 
particular  site  fj  is  occupied  is  either  one  or  zero;  nj  =  0,1.  The  average  occupation 
probability  <nj>  is  given  by  the  average  of  all  possible  configurations  of  the  particles 
each  weighted  by  its  statistical  probability  [14].  As  the  degree  of  order  changes  near 
a  phase  transition,  <nj>  also  changes. 

The  order  parameter  is  a  function  (or  set  of  functions)  that  describes  the 
average  degree  of  order  at  each  site  on  the  lattice  as  a  function  of  temperature.  In 
diffraction,  the  4>j  are  most  conveniently  defined  in  terms  of  the  Fourier  coefficients 
of  the  average  site  occupation  probability  <nj>.  The  number  of  components  j  of  the  the 
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order  parameter  is  equal  to  the  number  of  terms  in  this  expansion.  Consider  the 
example  in  Figure  3.  Figure  3a  shows  a  (2x1)  structure  on  a  square  lattice  which  has 
a  two  component  (vector)  order  parameter.  In  this  case  there  are  two  symmetry 
equivalent  orientations  of  the  (2x1 )  structure.  The  order  parameters  associated  with 
these  possibilities  are 


and 


i  v*  ikrri 

iv*  ik2ri 

*2 <Hj>  e 


(la) 


(lb) 


where  N  is  the  number  of  lattice  sites,  kj  *  aj*/2  and  k2  -  a2*/2  are  reciprocal 
lattice  vectors  characteristic  of  each  possible  (2x1)  lattice  as  shown  in  Figure  3b 
and  a  j*  and  a2*  are  the  reciprocal  lattice  unit  mesh  vectors  of  the  substrate.  For  a 
perfect  (2x1)  phase  4>|  =  1/2.  A  c(2x2)  phase,  Figure  2b,  on  this  same  substrate 
lattice  is  described  by  the  reciprocal  vector  q  and  has  a  single  component  (scalar) 
order  parameter 


*3=1  X<n|>e 


iqrj 


(1c) 


The  pair  correlation  function  <n(nj>  provides  more  detailed  information  about  the 
ordering  than  does  <t>  since  <njnj>  gives  the  probability  that  the  lattice  sites  at  rj  and 
r*j  are  simultaneously  occupied.  Since  n,  can  be  written  as  a  sum  of  its  average  value 
<nj>  and  its  instantaneous  fluctuation  6nj  from  that  average,  one  obtains 


<n(  nj>  =  <nj> 


2  + 


<SnjSnj>. 


(2) 


The  functional  form  of  <6nj6nj>  is  such  that  it  decreases  as  rj  and  rj  become  farther 
apart.  The  length  which  characterizes  this  fall-off  is  called  the  correlation  length  £ 
and  is  a  measure  of  the  distance  over  which  fluctuations  from  the  average  coverage 
are  <nj>  correlated.  Higher  order  correlation  functions  can  also  be  defined  but  little  is 
known  about  their  properties  for  real  systems. 
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2.2  Diffraction  from  a  Disordered  Lattice 

If  a  perfect  two-dimensional  lattice  is  illuminated  by  radiation  with  wavevector 
k0  (|k0l  *  2nA  ,  where  ^  is  the  wavelength),  diffracted  beams  are  observed  only  at 
scattered  wave  vectors  k  which  satisfy  the  Laue  conditions  [11,13],  If  disorder  is 
introduced  onto  the  lattice,  intensity  is  removed  from  the  beams  and  redistributed 
throughout  space.  In  this  case,  the  averaged  intensity  scattered  from  the  lattice  is 

I  =  lf(k,k0)|2  [lLR0  ♦  x(S)}  (3) 

where  |f(k,k0)l2  is  the  differential  scattering  cross  section  and  S  *  k  -  k0  is  the 
scattering  vector  [13].  Only  Sy,  kQj|  and  ky,  the  components  of  S,  k0  and  k,  which  are 
parallel  to  the  plane  of  the  lattice  are  important.  Equation  (3)  neglects  any  multiple 
scattering;  these  contributions  can  be  included  by  adding  terms  which  include  higher- 
order  correlation  functions  [15].  iLR0,  the  first  term  in  (3),  is  just  the  scattered 
intensity  that  would  be  obtained  for  a  perfect  lattice  with  a  particle  at  every  lattice 
site  but  reduced  by  the  factor  $2;  i.e. 

•lro  s  $2  6(s||  •  Ghk> 

where  are  the  reciprocal  lattice  vectors  of  the  two  dimensional  lattice. 

X(S),  the  second  term  in  (3)  is  due  to  fluctuations: 


In  analogy  with  magnetic  systems,  x(5)  is  called  the  susceptibility.  x(5)  is  a  peaked 
function  and  its  half-width  ASj|(  1  /2)  is  inversely  related  to  the  correlation  length  4 
by 

|AS||(  l/2)|  =  2ir/$.  (6) 

For  the  case  of  continuous  phase  transitions,  the  susceptibility  is  often  approximated 
as  a  Lorentzian  function  [  1 6, 1 7], 


X(S)  =  X  <«nj6nj 
U 
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2.3  Critical  Exponents  and  Universality  j 

j 

The  modern  theory  of  continuous  phase  transitions  postulates  that  certain  j 

entities  called  critical  exponents  provide  a  meaningful  description  of  a  system  near  a  | 

critical  point.  For  example,  if  F(T)  is  some  physical  quantity  of  interest  such  as  the 
order  parameter,  its  variation  near  the  critical  temperature  can  be  described  by 

F(t)  “  Foltl03  +  ...  (7) 


where  a)  is  the  critical  exponent,  t  =  (T  -  Tc)/Tc  is  the  reduced  temperature  and  F0  is 
an  amplitude.  Higher-order  terms  in  (7)  are  unimportant  as  long  as  t  is  near  enough  to 
zero.  Originally  it  was  hoped  that  the  critical  exponents  would  have  a  single  set  of 
values  [17],  However,  in  reality  there  seem  to  be  several  distinct  sets  which 
apparently  depend  upon  the  dimensionality  of  the  system  and  the  symmetry  of  the 
order  parameter  [  1 7],  Each  set  is  called  a  Universality  Class.  The  important 
universality  classes  for  two-dimensional  systems  are  listed  in  Table  1.  The  critical 
exponents  which  can  be  determined  by  a  diffraction  experiment  are  also  listed.  Table 
2  contains  a  compilation  of  critical  exponents  which  have  been  extracted  from 
experimental  data  from  reconstructed  clean  surfaces  and  chemisorbed  layers. 

The  temperature  dependence  of  the  order  parameter  is  described  by  the  critical 
exponent  fi; 


4>(t)  ^  4>0  ItlB  ♦ ...  T  <  Tc 


(8) 


=  0  T  >  T^ 

The  exponents  y  and  v  describe  the  behavior  of  x(AS||.t),  equation  (5).  The  amplitude 
X(0,t)  obeys 


X(0,t)  ^  Xo  itrv  ♦ ... 


(9) 


and  the  correlation  length  is  given  by 
4(t)  ^  £0ltrv  +  ... 


(10) 


where  £0  is  of  the  order  of  a  lattice  spacing. 
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Figure  4  is  a  schematic  illustration  of  the  variation  of  lLR0  and  x(5|j,t)  near  Tc 
predicted  by  equations  (8)  to  (10).  ILR0 falls  continuously  to  zero  at  Tc  while  x(O.t) 
has  a  singularity  at  Tc  and  decreases  smoothly  on  both  sides.  Far  below  Tc  only  the 
8-function  lLR0  contributes  to  the  line  profile.  As  T  increases  toward  Tc,  lLR0 
decreases  and  x(5||,t)  becomes  significant.  At  T  <<  Tc,  x(S||,t)  is  low  and  broad  but  it 
increases  in  magnitude  and  sharpens  until  at  Tc  it  is  the  only  contribution  to  the  line 
profile.  As  T  increases  above  TCj  x(S||,t)  again  decreases  and  broadens.  Thus,  an 
accurate  measurement  of  the  diffraction  line  profile  as  a  function  of  temperature  near 
Tc  can  be  analyzed  to  obtain  J3,  y  and  v.  This  is  the  theoretical  basis  of  most  studies 
of  second-order  phase  transitions. 

The  critical  exponent  a  is  traditionally  associated  with  the  specific  heat  anomaly 

C  =  C0ltl~a  ♦  ... .  (11) 

As  recently  pointed  out  by  Roelofs,  Einstein  and  Bartelt  [15],  a  also  describes  the 
behavior  of  the  “integrated"  diffraction  intensity  according  to  the  relation 

lint(t)  =  A0  +  A^iBjtl1**  +...  (12) 

where  I  jnt(t)  is  the  integral  of  the  intensity  profile  over  a  small  range  of  aSjj 
values  near  the  center  of  the  Brillouin  zone,  i.e.  at  Sjj  =  G^.  This  result  is  valid  as 
long  as  t  is  close  enough  to  zero  so  that  ASjj  is  larger  than  the  width  of  the 
susceptibility  a5||(1/2).  Equation  (12)  relies  on  the  fact  that  a  poor  LEED  instrument 
is  sensitive  predominantly  to  short-range  order.  This  observation  is  extremely 
important  because  it  opens  the  quantitative  study  of  critical  behavior  to  researchers 
with  standard  commercial  instruments.  Furthermore,  equation  (12)  is  valid  even  when 
multiple  scattering,  which  is  neglected  in  equation  (3),  is  important. 

2.4.  Critieria  for  a  Continuous  Order-Disorder  Phase  Transition:  Landau  Theory 

By  1937,  Landau  had  recognized  the  role  of  symmetry  in  continuous  order- 
disorder  phase  transitions  [18],  In  Landau's  theory,  the  thermodynamic  potentials 
$(9,T,4>)  are  expanded  in  powers  of  the  order  parameter  4>(0,T)  in  the  vicinity  of  the 
critical  point: 

4>(9,T,4>)  =  $0(e,T)  ♦  a(e,T)4>  +  A(e,T)4>2  ♦  B(e,T)4>3*  D(e,T)<J>4  v..  (1 3) 
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where  a(9,T),  A(0,T),  B(8,T)  and  O(0,T)  are  temperature  and  coverage  dependent  ; 

coefficients  which  can  themselves  be  expanded  in  powers  of  t  near  the  critical  point.  I 

The  order  parameter  4>  is  also  a  function  of  0  and  T  and  is  determined  by  minimizing  1 

4>(0,T)  at  particular  8  and  T.  Landau  and  Lifshitz  [  1 9]  determined  the  properties  of  the 
coefficients  of  equation  (13)  by  applying  a  set  of  four  symmetry  rules  which  are  now 
known  as  the  Landau-Lifshitz  Rules  [20]  and  are  discussed  in  detail  in  References  19  - 
22. 

The  Landau-Lifshitz  rules  make  it  possible  to  find  all  possible  second-order 
phase  transitions  (except  commensurate-incommensurate  transitions)  without 
knowledge  of  the  detailed  mechanism  of  the  ordering  process.  The  rules  however  do 
not  prove  that  a  given  system  will  actually  have  a  second-order  transition,  only  that 
such  a  transition  is  possible. 

Landau  theory  has  been  applied  to  two  dimensional  systems  and  a  complete 
classification  of  all  possible  order-disorder  phase  transitions  has  been  made  [21,22]. 

These  are  summarized  in  Table  1. 

The  procedure  for  predicting  the  universality  class  for  a  particular  system  is 
straightforward.  First,  the  components  of  the  order  parameter  must  be  identified. 

The  Hamiltonian  of  the  system  is  then  expanded  in  terms  of  these  components  and  the 
Landau-Lifshitz  rules  applied.  The  resulting,  symmetrized  Hamiltonian  is  called  the 
Landau-Ginzburg-Wilson  (LGW)  Hamiltonian.  This  LGW  Hamiltonian  is  then  comparied 
with  the  LGW  Hamiltonians  of  the  known  universality  classes  [20-22]  in  order  to 
identify  the  correct  universality  class.  Examples  are  given  in  references  20-23. 

Landau  theory  cannot  be  applied  too  near  the  critical  line  (0C,TC)  since  the 
expansion  of  $(0,T,4»)  in  equation  ( 13)  diverges  there.  However,  as  long  as  the 
root-mean-square  fluctuations  of  <J>  averaged  over  the  area  are  small  compared  to 
its  average  value,  the  arguments  given  above  are  expected  to  hold.  Ipatova  and  Kitaev 
[20]  have  pointed  out  that  this  condition  is  more  easily  met  in  diperiodic  three 
dimensional  systems  than  in  two  dimensional  systems.  They  suggest  that  diperiodic 
systems  may  provide  an  alternative  approach  for  real  surfaces. 

2.5  Finite  Size  Effects. 


The  characterization  of  continuous  phase  transitions  outlined  above  is  only 
precise  when  the  substrate  lattice  is  perfect,  e  g.  semi-infinite  in  extent  and  free  of 
imperfections.  At  best,  this  requirement  can  only  be  approximately  satisfied  by  real 
surfaces  since  imperfections  due  to  substrate  mosaic  structure,  steps,  dislocations 
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and  impurities  are  always  present.  The  influence  of  imperfections  on  the  critical 
properties  are  called  finite  size  effects. 

Usually,  finite  size  effects  become  important  whenever  the  temperature  is  close 
enough  to  Tc  that  the  correlation  length  £  becomes  comparable  to  the  spacing  between 
imperfections.  At  this  point  the  imperfections  begin  to  influence  the  fluctuations  and 
alter  the  critical  exponents  in  equations  (8H12).  The  effective  critical  temperature 
is  also  altered  [24,25],  The  surfaces  of  real  single  crystals  are  usually  only  free  of 
imperfections  on  the  scale  of  several  tens  of  nanometers  and  finite  size  effects  are 
important.  Quantitative  descriptions  of  the  way  in  which  finite  size  effects  influence 
critical  behavior  have  only  been  developed  for  a  few  special  geometries  and 
considerable  theoretical  work  is  still  needed.  Kleban  [23]  has  recently  reviewed  the 
status  of  this  area. 

3.  EXPERIMENTAL  METHODS 

Experimental  characterization  of  an  order-disorder  phase  transition  in  a  real 
surface  system  is  not  simple  for  several  reasons.  Sample  preparation  is  difficult  and 
not  wel  1-developed.  Accurate  absolute  measurement  and  control  of  substrate 
temperature  and  adsorbate  coverages  is  very  difficult.  The  time  available  for 
measurements  is  often  severely  limited  by  the  need  to  avoid  surface  contamination 
due  to  adsorption  of  residual  gases.  The  response  of  the  instrument  distorts  the 
measured  line  profiles  and  limits  the  sensitivity  to  long-range-order.  Other  processes 
also  contribute  to  changes  in  the  line  profile  and  are  not  easily  accounted  for  in  the 
analysis.  These  include  scattering  from  impurities,  extended  surface  defects,  the  bulk 
mosaic  structure  and  thermal  vibrations.  Furthermore,  in  the  case  of  LEED,  multiple 
scattering  is  a  major  contributor  to  the  total  intensity. 

Despite  this  long  list  of  experimental  concerns,  numerous  successful 
measurements  have  been  completed  and  substantial  progress  is  being  made  to 
overcome  most  of  these  outstanding  problems. 

3.1  Sample  Preparation  and  Surface  Imperfections. 

The  methods  currently  used  to  prepare  oriented,  single  crystal  surfaces  are  not 
highly  developed.  They  are  based  on  standard  metailographic  practices  [26]  and 
typically  include  mechanical,  chemical  and/or  electrochemical  polishing  sequences. 
Crystallographic  orientation  is  usually  determined  with  the  aid  of  X-ray  Laue  back 
reflection  photographs.  Since  this  method  is  only  accurate  to  3bout  ♦  0.2’,  a  typical 
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surface  has  atomic  height  steps  separated  on  average  by  not  more  than  a  few  tens  of 
nanometers.  Surfaces  of  this  quality  are  only  marginally  suited  for  studies  of 
long-range-order.  More  accurate  X-ray  orientation  procedures  requiring  measurement 
of  line  profiles  and  relative  intensities  can  significantly  reduce  the  degree  of 
misorientation  to  less  than  0.05’  [27]  and  need  to  be  more  widely  employed.  Surface 
polishing  procedures  also  need  to  be  more  carefully  developed  to  minimize  degree  of 
surface  disorder. 

The  surface  cleaning  procedures  employed  inside  the  ultra-high  vacuum  analysis 
chamber  also  have  considerable  potential  for  introducing  surface  imperfections.  The 
"standard"  method  is  sputter  cleaning  by  bombardment  with  energetic  inert  gas  ions. 
Sputtering  is  a  violent  process  and  creates  extensive  damage  usually  in  the  form  of 
randomly  distributed  step  arrays  but  sometimes  also  creates  irreversible  macroscopic 
changes  including  formation  of  cones  [28],  The  degree  to  which  vacuum  annealing 
restores  a  defect  free  surface  is  not  well  understood.  Chemical  cleaning  processes 
[29]  involving  reaction  with  gases  are  intuitively  appealing  but  no  studies  of  their 
effects  on  the  surface  long-range-  order  have  been  made.  Unfortunately,  chemical 
reactions  may  not  be  able  to  remove  all  impurities. 

Surface  spectroscopy  techniques  such  as  Auger  electron  spectroscopy  [13,30] 
seldom  have  practical  sensitivities  better  than  about  0.1%  of  a  monolayer.  Thus,  the 
mean  separation  of  surface  impurities  can  easily  be  smaller  than  10  nanometers. 
Adsorption  of  impurities  from  the  residual  gas  is  also  an  important  source  of 
contamination  for  many  materials.  For  example,  for  a  pressure  of  ICT12  Torr  and 
unity  stickinq  probability,  less  than  one-half  hour  is  required  to  accumulate  about 
0.1%  of  a  monolayer  of  impurity.  This  corresponds  to  an  impurity  atom  separation  of 
about  7-8  nm. 

The  extent  to  which  uncontrolled  defects  and  impurities  influence  studies  of 
critical  phenomena  by  introducing  finite  size  effects  is  not  known  at  present.  The 
improved  diffraction  equipment  needed  to  better  characterize  the  starting  surface  is 
only  just  becoming  available  [31]. 

Studies  of  critical  behavior  of  chemisorbed  layers  on  surfaces  have  relied  on 
Auger  electron  spectroscopy  (AES)  to  determine  the  relative  surface  coverage  of  the 
adsorbate.  For  submonolayer  coverages,  the  AES  signal  is  directly  proportional  to  the 
coverage  as  long  as  the  electronic  structure  of  the  adsorbate  does  not  change 
drastically  with  coverage.  Thus,  the  relative  coverage  can  be  determined  with  high 
accuracy.  However,  AES  cannot  measure  the  absolute  coverage  to  better  than  twenty 
or  thirty  percent  so  that  another  method  is  always  needed  to  calibrate  the  coverage 
scale.  This  method  has  been  LEED  in  all  of  the  experimental  studies  of  overlayer 
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critical  phenomena  discussed  in  the  next  section.  For  example,  the  coverage  scale  in 
Figure  1  was  set  by  assuming  that  the  half-order  diffraction  beams  of  the  p(2  x  2) 
phase  reach  maximum  intensity  at  precisely  one-quarter  monolayer  [5].  This  method 
can  easily  be  in  error  by  five  to  ten  percent.  Kinetic  effects  may  limit  the  removal  of 
domain  boundaries.  Adsorption  near  impurities  and  defects  may  have  different 
properties  than  on  the  perfect  surface  [23,32],  Even  absolute  methods  such  as  nuclear 
reaction  analysis  [33]  and  x-ray  fluoresence  [34]  are  not  more  accurate  than  a  few 
percent.  Fortunately,  for  many  properties  of  interest,  it  is  sufficient  to  know  only 
relative  coverages. 

Imperfections  give  rise  to  finite  size  effects  which  limit  the  smallest  value  of  t 
for  which  the  power  law  behavior  given  in  equations  (8)  -  (12)  will  be  valid.  The 
smallest  value  of  t  attained  so  far  is  about  0.005  [35].  This  has  the  consequence  that 
only  very  small  ranges  of  t  have  been  used  to  test  equations  (8)  -  (12),  typically  about 
one  order  of  magnitude.  Thus,  at  present,  surface  imperfections  provide  one  of  the 
most  significant  limitations  to  our  ability  to  test  universality  in  model 
two-dimensional  systems. 

3.2  Diffraction  Instrumentation 

Low-Energy  Electron  Diffraction  (LEED)  is  the  only  technique  that  has  been  used 
to  study  critical  phenomena  at  continuous  order-disorder  transitions  for  systems  of 
interest  in  this  review.  Atom  diffraction  [36]  and  X-ray  Diffraction  [37]  would  also  be 
suitable  but  have  not  yet  been  employed,  in  LEED,  not  only  do  sub-surface  atomic 
layers  contribute  to  the  signal  but  multiple  scattering  effects  may  be  too  large  to  be 
ignored  [38,39],  This  makes  complete  analysis  of  LEED  intensities  complex. 

An  ideal  instrument  would  illuminate  the  sample  with  a  plane  wave  of 
monochromatic  radiation  and  the  detector  would  accept  only  elastically  scattered 
radiation  with  a  specific  wavevector  k.  Real  instruments  do  not  meet  these 
requirements  [40],  Deviations  from  ideal  behavior  can  be  described  by  the  instrument 
response  function  T(5)  [41],  If  1(5)  represents  the  diffracted  intensity  in  the  absence 
of  instrumental  effects,  then  the  measured  intensity  J(5),  is  given  by 

J(5)  =  T(5)  *  1(5)  (14) 

where  *  represents  a  convolution  product.  In  principle,  most  features  of  T(5)  can  be 
determined  experimentally  ]f  a  sufficiently  perfect  crystal  is  available  and  the 
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determination  is  carried  out  at  a  temperature  which  is  low  enough  that  thermal 
diffuse  scattering  due  to  phonons  is  negligible. 

T(S)  is  often  characterized  by  its  Fourier  transform  t(R);  t(R)  is  peaked  at  R  =  0 
and  falls  off  to  zero  at  large  R.  The  transfer  width  R,/2  is  the  full-width  at  half¬ 
maximum  of  t(R)  and  is  used  as  a  measure  of  the  degree  to  which  a  particular 
diffractometer  is  sensitive  to  long-range-order.  For  an  ideal  instrument,  R1/2  would 
be  infinite  but  for  most  LEED  instruments  currently  in  use,  R1/2  is  small,  typically  in 
the  range  5-20  nm.  Since,  correlation  lengths  can  be  several  hundred  Angstroms  in 
two-dimensional  systems  [6],  a  small  value  of  Rl/2  can  make  it  difficult  to 
distinguish  certain  types  of  fluid  from  solids.  Recently,  several  new  LEED 
diffractometer  designs  have  been  realized  and  provide  substantially  improved 
instrumental  resolution  [31,40]. 

Whenever  the  correlation  length  i  is  comparable  to  R1/2,  it  is  essential  that  the 
effect  of  T(S)  on  the  measured  line  profile  be  included  in  the  analysis.  Ideally,  this 
would  be  done  by  direct  deconvolution  of  equation  (14)  to  extract  KS),  however  this  is 
difficult  in  practice  both  because  of  uncertainties  in  the  measured  T(S)  and  because  of 
the  finite  range  of  5  over  which  it  can  be  determined.  Instead,  the  normal  procedure  is 
to  assume  that  the  width  of  a  profile  measured  far  below  Tc  is  entirely  due  to  the 
instrument  response.  In  the  vicinity  of  Tc,  the  measured  profile  may  then  be  fit  by 
assuming  a  particular  form  for  the  fluctuation  contribution  x(S||,t),  equation  (5), 
usually  x(S||,t)  is  taken  to  be  Lorentzian  [16],  The  amplitude  Xo  and  width  A5||(  1/2)  of 
the  Lorenztian  plus  the  intensity  of  the  low-temperature  line  profile  are  used  as 
adjustable  fitting  parameters.  The  temperature  dependence  of  each  of  these 
parameters  is  then  fit  to  equations  (8)-(  12),  to  determine  the  critical  exponents.  A 
constant  term  is  usually  included  in  the  analysis  to  approximately  account  for  a 
"background". 

3.3.  Additional  Considerations. 

There  are  several  additional  factors  which  can  change  the  shape  or  temperature 
dependence  of  a  line  profile.  For  example,  the  experimental  configurations  normally 
employed  do  not  measure  the  diffraction  line  profile  along  S||,  but  along  some  curved 
trajectory  in  reciprocal  space  [42],  This  can  result  in  an  underestimate  of  A5||(  1  /2) 
by  as  much  as  fifty  percent  and  thus  in  an  overestimate  of  the  correlation  length  £. 

In  addition,  the  equations  presented  in  Section  2  do  not  incorporate  the  effects  of 
all  physical  processes  which  contribute  to  the  diffracted  intensity  and  its  line  profile. 
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One  of  these  is  the  temperature  dependence  due  to  the  Debye-Waller  factor  e-21*1  which 
is  a  function  of  both  5  and  T; 

2M  =*  <(S-u)2> 

where  u  is  the  root-mean-square  vibrational  amplitude  of  the  surface  atoms  and  is 
proportional  to  T  in  the  high  temperature  limit.  2M  can  be  determined  experimentally 
from  the  slope  of  In  I  vs.  T  measured  for  T  far  away  from  Tc.  It  is  typically  in  the 
range  0.5-3.0  for  most  materials  and  is  a  substantially  larger  correction  for  surfaces 
than  for  the  bulk  not  only  because  of  the  larger  vibrational  amplitudes  of  surface 
atoms  but  also  because  of  the  absense  of  translational  symmetry  normal  to  the 
surface  [43,44], 

The  intensity  removed  from  a  diffraction  beam  by  the  Debye-Waller  factor  is 
redistributed  throughout  the  Brillouin  zone  as  one-  and  multi-phonon  thermal  diffuse 
scattering  (TDS)  [43,44],  The  intensity  of  the  one-phonon  scattering  has  its  maximum 
near  the  Bragg  peak.  This  additional  peaked  intensity  can  cause  the  measured  line 
width  to  be  increased  by  twenty  to  thirty  percent  [42]!  The  effects  of  TDS  increase  at 
higher  T  and  for  large  values  of  S.  No  experimental  analysis  of  critical  behavior  has 
yet  attempted  to  include  the  effects  of  TDS. 

Multiple  scattering  is  very  strong  for  low  energy  electrons  and  diffracted 
intensities  cannot  be  calculated  using  the  single  scattering  theory  as  assumed  in  the 
derivation  leading  to  equation  (3).  Consequently,  one  must  ask  whether  multiple 
scattering  can  alter  the  line  profiles  sufficiently  to  invalidate  the  use  of  equation  (3) 
for  the  analysis  of  real  LEED  data.  Unfortunately,  we  do  not  have  clear  answer  to  this 
question  at  the  present  time.  At  least  some  experiments  have  attempted  to  use 
diffraction  conditions  which  might  minimize  the  contributions  of  multiple  scattering. 
However,  the  fact  that  critical  phenomena  are  studied  using  non-specular  beams 
means  that  at  least  two  Bloch  waves  can  be  excited  and  must  be  treated  self- 
consistently  in  any  analysis  [45],  Recent  results  of  conformal  mapping  analysis 
indicate  that  multiple  scattering  might  even  dominate  the  line  profiles  at  Tc  [46].  A 
recent  theoretical  analysis  by  Moritz  and  Lagally  [39]  suggests  that  multiple 
scattering  cannot  be  ignored.  However,  Bennett  and  Webb  [47]  have  made  line  profile 
measurements  on  several  beams  and  they  found  consistent  results.  Since  the  relative 
contributions  from  multiple  and  single  scattering  will  vary  strongly  from  diffraction 
beam  to  diffraction  beam,  this  result  gives  some  support  for  neglect  of  multiple 
scattering  in  line  profile  analysis.  Comparison  of  profiles  measured  for  the  same 
system  with  both  LEED  and  with  X-ray  diffraction  could  directly  test  for  the 
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importance  of  multiple  scattering  because  X-ray  intensities  are  described  by  the 
single  scattering  theory.  Such  measurements  are  planned  [48], 

Finally,  it  should  also  be  pointed  out  that  the  method  of  Bartelt,  Roelofs  and 
Einstein  [15]  for  determination  of  the  specific  heat  critical  exponent  a  does  not  seem 
to  be  sensitive  to  multiple  scattering  effects.  This  fact  makes  this  method 
additionally  attractive  to  the  experimentalist. 

4.  EXPERIMENTAL  STUDIE5  OF  CRITICAL  BEHAVIOR. 

Davisson  and  Germer  probably  made  the  first  observation  of  a  surface  order- 
disorder  transition  in  their  classic  paper  on  electron  diffraction  [49].  The  first 
modern  observation  was  by  McRae  for  0  on  Ni(  1 00)  [50];  this  system  was  later  shown 
to  be  complicated  by  the  onset  of  0  dissolution  into  the  bulk  for  T  ~  Tc  [51].  Perhaps 
the  first  experiment  designed  specifically  to  study  a  surface  order-  disorder 
transition  was  a  LEED  study  made  on  Cu3Au(  100)  by  Sundaram  et  al.  [52]. 

Only  a  few  careful  studies  of  critical  behavior  on  clean  and  chemisorbate  covered 
surfaces  have  now  been  carried  out.  These  are  reviewed  in  this  section. 

4. 1 .  Reconstructions  on  Clean  Surfaces. 

4.1.1.  Au(  1 1 0)  ( 1  x2):  Realization  of  the  tsina  Universality  Class. 

The  first  attempt  to  extract  a  critical  exponent  for  a  surface  system  was  made 
by  Wolf  et  al.  [53]  for  the  Au(  1 1 0)  ( 1  x2)  order-disorder  transition.  Although  they 
corrected  the  data  for  the  instrument  response  function,  they  did  not  include  critical 
fluctuations  in  the  analysis.  This  resulted  in  a  larger  value  for  13  than  currently 
accepted  for  this  system. 

A  complete  analysis  of  this  transition  was  made  by  Campuzano  et  al.  [16]  and 
confirmed  the  earlier  prediction  by  Bak  [54]  that  the  transition  should  be  in  the  simple 
Ising  universality  class.  The  line  profiles  of  the  (10),  (l, -1/2)  and  (1,7)  beams  were 
measured  along  the  <01>  direction  as  a  function  of  temperature  for  normally  incident, 

20  eV  electrons.  The  line  profile  measured  in  the  fully  ordered  (1x2)  phase  for  T  <  Tc 
is  a  measure  of  lLR0,  equation  (4),  and  could  be  fit  by  a  Gaussian  plus  a  uniform 
background  term.  The  width  of  this  Gaussian  was  substantially  larger  than  that  of  the 
instrument  response  function  and  indicated  that  finite  size  effects  limited  the  extent 
of  the  ( 1  x2)  surface  long-range  order  to  less  than  1 5  nm.  After  correction  for  the 
measured  Debye-Waller  factor,  the  experimental  line  profiles  were  fit  at  each 
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temperature  by  assuming  them  to  be  a  linear  combination  of  lLR0  and  the  susceptibility 
X(S||,t),  equation  (5),  plus  a  uniform  background.  The  critical  exponents  3,  y  and  v 
were  extracted  from  the  temperature  dependence  of  1LR0>  x0  and  ASj|  respectively 
assuming  the  forms  in  equations  (8)-(12).  The  resulting  fits  are  shown  in  Figure  (5a) 
and  the  exponents  are  tabulated  in  Table  2.  Tc  was  found  to  be  650  ♦  7  K. 

The  consequences  of  finite  size  effects  due  to  the  limited  perfection  of  the  (1x2) 
phase  are  reflected  in  both  lLR0(T)  and  Xo^  as  shown  in  Figure  5b.  ILR0  decreases 
slightly  less  steeply  than  predicted  by  the  classical  Onsager  model  [55],  This  causes 
the  exponent  13  to  be  slightly  larger  than  the  predicted  value.  Xo^  is  clearly  peaked 
above  Tc  in  a  way  similar  to  that  predicted  by  Fisher  and  Burford  [56],  As  shown  in 
Table  2,  the  measured  exponents  are  in  good  agreement  with  the  predicted  values 
except  for  13  which  is  slightly  too  large. 

Clark  et  al.  (35]  have  measured  the  exponent  of  the  specific  heat  anomaly  to  be  a 
=  0  i  0.05.  They  used  the  method  proposed  by  Bartelt,  Einstein  and  Roelofs  [  1 5]  which 
exploits  the  fact  that  the  integrated  intensity  measured  by  a  poor  LEtD  instrument  has 
lost  its  sensitivity  to  long-range  order  and  is  described  by  equation  (12).  They  found 
the  critical  temperature  to  be  690  ♦  7  K,  a  value  that  is  higher  than  that  reported 
Campuzano  et  al.  [16].  This  difference  is  due  to  the  substantially  better  surface  order 
of  the  crystal  used  by  Clark  et  al.  and  is  in  good  agreement  with  theoretical 
predictions  for  the  Ising  model  [24,25]. 

Because  Au(l  10)  is  inert  to  the  background  gases  normally  found  in  ultra-high 
vacuum  systems,  very  long  times  are  available  for  measurements.  Thus,  Au  is 
expected  to  be  a  good  model  system  for  more  extensive  studies  of  continuous  order- 
disorder  transitions  in  two  dimensions.  Finite  size  effects  can  be  introduced  by  using 
stepped  surfaces  or  intentionally  adding  impurities.  The  role  of  multiple  scattering 
and  other  corrections  to  the  data  such  as  TD5  can  be  explored  in  detail.  Studies  of  the 
kinetics  of  ordering  will  also  be  possible  since  the  time  scale  over  which  the  ordering 
takes  place  is  reasonably  slow. 

The  clean  lr(  1 1 0)  and  Pt(  1 1 0)  surfaces  also  have  a  ( 1  x2)  reconstruction  [57]  and 
are  thus  expected  to  have  order-disorder  transitions  which  belong  to  the  Ising 
universality  class.  Experiments  to  confirm  this  prediction  remain  to  be  carried  out. 

4.1.2.  W(001  )c(2x2):  The  XY  Model  with  Cubic  Anisotropy. 

The  clean  W(00 1 )  surface  reconstructs  to  a  c(2x2)  phase  with  a  critical 
temperature  of  21 1  K  [57-59],  The  diffraction  pattern  of  the  c(2x2)  phase  has  p2mg 
symmetry  due  to  small  displacive  translations  of  W  surface  atoms  along  <1 1>  and  <1  T> 
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directions.  Since  the  unreconstructed  ( 1  x  1 )  surface  has  p4mm  symmetry,  the  only 
irreducible  representation  which  can  have  a  continuous  order-disorder  transition  has  a 
three  component  order  parameter  with  the  wavevectors  q,  k  j  and  k2  shown  in  Figure 
3.  The  order  parameter  $3  corresponds  to  the  c(2x2)  phase  [60]  and,  because  (2x1 )  and 
(1x2)  phases  do  not  occur,  $|  =  <$>2  =  0.  The  possibility  of  two  independent  sets  of 
atomic  displacements  along  <  1 1  >  or  <  1  T>  directions  means  that  $3  has  two 
components 


and 


4>a  =  AacosqaT 
4>b  =  Abcosqb-r. 


The  corresponding  LGW  Hamiltonian  has  the  form 

HLqw  ^  (r/2)(4>a2  +  4>b2)  +  u(4>a2  +  4>b2)2  +  w(4»a4>b)^- 

The  third  term  breaks  the  isotropy  of  the  first  two  and  restricts  the  displacements  to 
lie  along  either  the  <  1 1  >  or  <  1  T>  directions.  This  restriction  to  orthogonal  directions 
is  refered  to  as  the  "cubic  anisotropy".  The  coupling  coefficient  w  must  be  positive  to 
insure  that  HL6W  has  a  minimum  value  when  the  displacemants  go  to  zero  at  the 
critical  temperature.  This  Hamiltonian  is  identical  to  that  of  the  XY  Model  with  cubic 
anisotropy  universality  class  [21,22],  This  class  is  interesting  because  its  critical 
exponents  are  non-universal,  that  is,  they  depend  upon  the  strength  of  the  interactions 
as  indicated  by  the  results  of  model  calculations  [6 1  ]  shown  in  Figure  6.  The  value  of 
P  is  predicted  to  be  very  near  the  two-dimensional  Ising  value  1/8  except  for  very 
small  w,  p  diverges  rapidly  to  infinity  at  w  =  0. 

Wendelken  and  Wang  [59]  have  measured  p  =  0. 1 4  ±  0.04,  which  is  within 
experimental  uncertainty  of  the  Ising  value  indicating  that  the  coupling  coefficient  w 
is  not  too  near  zero.  They  have  also  measured  p  on  a  surface  oriented  3.25*  off  the 
(001)  surface  to  produce  3.0  nm  wide  terraces  with  step  edges  parallel  to  [010]. 

These  finite  size  effects  caused  the  critical  temperature  to  incresase  to  217  K  and  p 
to  decrease  to  0.05  ±0.01.  This  large  change  in  p  illustrates  the  importance  of  finite 
size  effects  in  surface  studies  and  emphasizes  the  need  for  a  well-characterized 
substrate. 

Mo(00 1 )  undergoes  a  similar  transition  [62]  but  its  critical  behavior  has  not  been 
studied  in  detail. 
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4.1.3.  Other  Clean  Surface  Reconstruction  Systems. 

Order-disorder  phase  transitions  have  been  observed  for  several  other  clean 
surfaces.  However,  these  transitions  are  either  believed  to  be  first  order  or  to  mirror 
an  underlying  bulk  transition. 

Bennett  and  Webb  [47]  have  shown  that  the  clean  Si(  111)  (7x7)  surface  undergoes  _ 
a  reversible  order-disorder  transition  with  0  =  0.11  ♦  0.02  and  Tc  ~  1 140  K.  However, 
a  consistent  interpretation  of  these  results  in  terms  of  a  continuous  transition  was 
not  possible  since  the  data  showed  no  indication  of  the  expected  broadenirg  of  the  line 
profiles  due  to  critical  fluctuations.  This  transition  is  not  expected  to  be  continuous. 
The  Si(l  1 1)  (1x1)  surface  has  symmetry  with  six  classes  of  operators.  Thus,  the 
order  parameter  of  the  (7x7)  phase  can  contain  at  most  six  components  if  the 
transition  is  to  be  continuous.  However,  as  pointed  out  by  Willis  [60],  the  order 
parameter  may  have  up  to  48  components;  i.e.  the  number  of  new  diffraction  beams  per 
( lxl )  substrate  beam.  Although  the  actual  number  of  components  may  be  smaller  than 
48,  it  seems  unlikely  that  the  Landau  rules  can  be  satisfied  for  Si(7x7). 

The  bulk  disordering  of  CU3AU  is  a  classic  example  of  a  discontinuous  transition. 
A  transition  also  takes  place  in  the  same  temperature  range  on  the  (001)  surface  but 
with  modif ied  form.  This  surface  transition  has  been  studied  several  times  since  the 
early  experiments  of  Sundaram  et  al.  [52]  yet,  other  than  the  fact  that  it  accompanies 
the  bulk  transition,  its  detailed  character  remains  uncertain.  Jamison  et  al.  [63]  and 
Alvarado  et  al.  [64]  provide  evidence  using  polarized  LEED  that  this  transition  is 
continuous  although  the  detailed  intensity  and  polarization  variations  they  report  are 
substantially  different.  On  the  other  hand,  the  data  of  McRae  and  Malic  [65]  point  to  a 
discontiuous  transition.  Theoretical  models  by  Kumar  and  Bennemann  [66]  indicate 
that  either  case  can  be  possible  depending  upon  the  details  of  the  interatomic 
potentials  near  the  surface.  Substantial  work  needs  to  be  done  to  clarify  the 
experimental  situation. 

The  critical  behavior  of  the  surface  magnetization  of  the  Ni(00 1 )  surface  is 
another  surface  phase  transition  whose  critical  behavior  follows  the  bulk  transition. 
The  critical  exponent  p  has  the  value  0.825  and  lies  between  the  three  dimensional 
Ising  and  XY  Model  values  [67], 

4.2.  Continuous  Transitions  in  Adsorbed  Overlavers. 

4.2.1.  0/Ni(  111)  p(2x2) 
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The  order-disorder  transition  of  the  0  p(2x2)  phase  on  Ni(  1 1 1 )  is  the  first 
example  of  critical  exponent  determination  for  a  chemisorbed  layer  [68-70],  Line 
profiles  of  the  (0,1/2)  beam  were  measured  with  an  external  photometer.  The  profiles 
were  analyzed  using  a  method  similar  to  that  described  above  for  Au(  1 10)  except  that 
the  effects  of  the  broader  instrument  response  function  were  included.  The  resulting 
exponents  are  shown  in  Table  2.  The  lattice  of  adsorption  sites  on  Ni(  111)  has 
triangular  symmetry  since  the  oxygen  atoms  are  known  to  occupy  the  threefold 
hollows  so  that,  at  first  glance,  Table  1,  this  transition  is  expected  to  be  in  the  so- 
called  four-state  Potts  universality  class  [21,22], 

The  measured  exponents  are  clearly  much  larger  than  the  predicted  values,  in 
fact,  within  the  experimental  uncertainties,  they  agree  with  the  Ising  exponents. 
Roelofs  et  al.  [68]  were  able  to  show  that  this  discrepancy  cannot  be  explained  by  the 
logariathmic  corrections  to  the  four-state  Potts  model  which  were  not  included  in 
their  original  analysis  of  the  data.  Several  explanations  for  this  have  been  suggested. 
Schick  [71]  pointed  out  that  there  are  two  types  of  adsorption  sites  of  the  (1 1 1) 
surface  of  a  face-centered-cubic  crystal— those  which  would  be  occupied  by  the  next 
layer  of  nickel  atoms  in  the  normal  abcabc  ...  stacking  sequence  and  those  sites  which 
would  be  hexagonal-close-packed-like  with  ababcabc  ... .  The  four-state  Potts  model 
results  only  if  there  is  a  significant  energy  difference  between  these  two  types  of 
sites.  If  there  is  no  difference,  the  sites  become  equivalent  forming  a  honeycomb 
lattice.  In  this  case,  the  transition  could  be  in  the  Heisenberg  model  with  cubic 
anisotropy  until  very  near  Tc  where  4-state  Potts  behavior  would  dominate.  The 
Heisenberg  model  with  cubic  anisotropy  has  variable  exponents  which  can  have  values 
very  near  those  of  the  Ising  model  [72], 

A  second  explanat.on  is  that  defects  in  the  form  of  impurities,  random  defects 
and  edge  effects  can  change  the  critical  behavior  from  that  expected  for  a  perfect 
infinite  system  [73,74],  For  example,  calculations  have  demonstrated  that  random 
impurities  dramatically  change  the  exponents  to  values  in  agreement  with  the 
measurements  [75],  The  nickel  substrates  used  in  these  experiments  were  more 
carefully  prepared  than  is  typically  the  case  in  surface  experiments  making  it  unlikely 
that  more  perfect  surfaces  can  be  easily  obtained  v’thout  using  more  sophisticated 
X-ray  orientation  methods.  Furthermore,  in  any  real  system  the  mosaic  structure  and 
equilibrium  concentration  of  vacancies  and  self-adatoms  can  be  substantial.  Thus,  if 
this  latter  explanation  is  true,  critical  exponents  extracted  from  surface  systems 
would  have  to  be  interpreted  with  great  caution. 

A  third  explanation  is  that  other  effects  such  as  multiple  scattering,  thermal 
diffuse  scattering  or  variations  in  the  effective  nickel  scattering  factors  are 
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important  [40],  Repeating  the  analysis  at  significantly  different  diffraction 
geometries  provides  one  way  to  test  for  these  effects. 

4.2.2.  Ao(1  10)  c(2x2)  Cl. 

Chlorine  occupies  the  fourfold  hollow  sites  on  the  square  lattice  of  the  Ag(lOO) 
surface  and  forms  an  ordered  c(2x2)  structure  [76].  The  order-disorder  transition 
should  thus  be  in  the  I  sing  Universality  Class.  A  recent  experimental  study  by  Tayler 
et  al.  [77]  shows  that  this  indeed  seems  to  be  the  case.  However,  only  one  critical 
quantity  was  extracted  from  the  data. 

The  only  critical  line  observed  for  this  system  is  that  separating  the  low 
coverage  disordered  phase  from  the  higher  coverage  c(2x2)  phase.  Since  this  line  is 
nearly  vertical,  Cl/Ag(  1 00)  is  an  example  of  a  system  in  which  a  critical  coverage,  0C, 
rather  than  a  critical  temperature,  must  be  measured.  As  discussed  by  Tayler  et  al., 
the  critical  behavior  of  the  intensity  is  decribed  by  an  effective  exponent  p/(l-a). 

The  value  of  p/(  1  -ot)  is  found  to  be  0. 1 2  +  0.03  in  good  agreement  with  the  Ising  value. 
Table  2.  The  critical  coverage  is  0C  =  3.94+  0.007  independent  of  temperature  up  to 
the  onset  of  Cl  desorption.  This  value  is  slightly  higher  than  the  hard  square  Ising 
mode]  value  of  0.368  probably  due  to  a  small  next  nearest  neighbor  repulsion. 

Although  the  expected  broadening  of  the  line  profiles  was  observed,  the  critical 
exponents  y  and  v  were  not  extracted  from  the  data.  5ince  there  is  some 
disagreement  about  the  behavior  of  Cl  on  Ag(  100)  [76],  determination  of  these 
additional  exponents  is  important  to  confirm  its  universality  classification. 

4-2.3.  W(1 10)  p(2x1)-H  and  d  (2x2)-H. 

Lyuksyutov  and  Fedorus  [78]  have  estimated  a  value  of  (3  from  their  measurements 
of  the  T-dependence  of  the  fractional  order  diffraction  beams  on  this  adsorption 
system.  They  do  not  report  line  profile  measurements  and  did  not  attempt  to  account 
for  the  effects  of  instrument  response  or  critical  fluctuations  in  their  analysis.  In 
order  to  account  for  possible  finite  size  effects  on  their  results,  they  assumed  a 
Gaussian  distribution  of  Tc  with  a  width  aT  in  their  fits  to  the  measured  data.  In 
total  three  parameters  (Tc,  13,  aT)  were  used  in  the  fits.  The  large  errors  in  the  values 
of  p,  Table  2,  are  due  to  the  strong  interdependence  of  I|_rq  on  P  and  Tc- 

The  bonding  site  of  H  on  W(  1 1 0)  is  the  subject  of  some  controversy  since  it  has 
only  recently  been  realized  that  the  substrate  reconstructs  [79],  If  it  is  assumed  that 
H  bonds  either  linearly  on  top  of  W  atoms  or  in  equivalent  long  bridge  sites,  the 
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symmetry  is  the  same  as  the  bulk,  c2mm,  and  these  phase  transitions  are  predicted  to 
belong  to  the  XY  Model  with  cubic  anisotropy  univesality  class.  As  shown  in  Figure  6,  p 
should  be  very  near  the  Ising  value  0.125  unless  the  cubic  interaction  w  is  very  small, 
p  diverges  rapidly  to  infinity  at  w  =  0.  Lyuksyutov  and  Fedorus’s  low  coverage  values 
are  in  agreement  with  this  model  but  the  critical  exponent  values  seem  too  high. 

Substantially  more  work  is  required  on  this  system.  In  particular,  a  more 
complete  experimental  study  of  the  critical  behaviors  essential.  Since  the  fractional 
order  beams  due  to  ordered  H  overlayers  are  weak  even  in  the  presence  of  substrate 
reconstruction,  such  studies  are  expected  to  be  difficult.  Measurements  should 
concentrate  near  the  so-called  Kosterlitz-Thouless  point  where  the  critical  lines 
separating  the  (2x1),  (2x2)  and  fluid  phases  come  together. 

4.2.4.  Oxygen  on  W(  1 12), 

Wang  and  Lu  [80]  have  studied  the  critical  behavior  of  the  p(2xl)  structure  formed 
by  oxygen  chemisorbed  on  the  WM  12)  surface.  The  critical  exponents  extracted  from 
the  experimental  data,  Table  2,  support  the  Ising  universality  classification  obtained 
from  symmetry  arguments.  Above  Tc,  measured  profiles  were  fit  in  the  way  described 
above  for  Au(  1 10)  [16]  to  obtain  y  andu  The  resulting  values,  y  =  1.79  t  0.14  and  v  = 

1 .09  i  0. 1 1 ,  are  in  good  agreement  with  the  Ising  values  of  1  /8  and  I ,  respectively. 

The  exponent  B  =  0.13  i  0.0 1  was  extracted  without  correcting  the  data  for  the 
instrument  response  using  a  method  proposed  by  Lyuksyutov  and  Fedorus  [78],  This 
method  cannot  be  applied  too  near  Tc  since  it  does  not  account  for  the  contributions 
due  to  fluctuations.  This  greatly  restricts  the  range  of  t,  Wang  and  Lu  found  0.016  <  t 
<  0.083.  A  reevaluation  of  their  data  using  the  BER  method  [15]  should  also  allow  the 
exponent  a  to  be  determined. 

As  is  the  case  for  most  experimental  studies  of  phase  transitions,  dat3  was 
reported  for  only  a  single  fractional  order  beam. 

4.2.5.  Group  1  a  and  2a  Metals  on  W(Q  1 1)  and  Mo(0 11 ). 


Fairly  extensive  studies  of  the  phase  diagrams  of  the  Group  la  and  2a  metals 
adsorbed  on  the  (Oil)  surfaces  of  W  and  Mo  have  been  carried  out  [  1 ,8 1  -  83],  Ba 
adsorbed  on  W(0 1 1 )  is  a  typical  example  [81].  A  series  of  coincidence  phases  is 
observed  at  low  coverages  and  undergo  first-order  transitions  upon  heating.  However, 
a  commensurate  (3x2)  phase  at  0  “  0.83  appears  to  disorder  continuously  upon  heating 
with  a  critical  temperature  near  1 25  K.  Fedorous  and  Gonchar  [8 1  ]  measured  the 


UNERTL/REVIEW 


intensities  of  fractional  order  beams  during  this  transition  using  a  detector  aperature 
approximately  equal  to  the  full-width  at  half-maximum  of  diffraction  beams  in  the 
well-ordered  low-temperature  (3x2)  phase.  They  erroneously  attempted  to  extract  the 
critical  exponent  13  from  this  data  by  ignoring  the  effect  of  their  poor  instrumental 
response  function.  An  appropriate  analysis  should  use  the  BER  method  [16).  We  have 
applied  such  an  analysis  to  the  data  in  Figure  3  of  reference  81  using  the  assumptions 
described  in  reference  55.  We  find  that  a  is  near  zero  for  the  (3x2)  order-disorder 
transition  indicating  that  it  must  belong  to  either  the  simple  Ising  or  XY  with  cubic 
anisotropy  universality  classes. 

4.2.6.  Other  Adsorbate  Systems. 

The  phase  diagrams  of  a  large  number  of  other  adsorbate  systems  have  been 
characterized.  A  few  of  these  seem  to  have  continuous  transitions  and  should  be  the 
subject  of  future  studies.  Some  examples  include:  Se/Ni(00 1 )  [5],  5/Ni(  111)  [84], 
Te/Ni(  111)  [35],  Bi/Cu(  111)  [86],  H/W(00 1 )  [60],  H/Mo(00 1 )  [60],  S/Mo(  1  1 0)  [87], 

4.2.7.  Two  Component  Adsorbed  Lavers 

Multi-component  adsorbed  layers  should  offer  even  richer  behavior  than  the  single 
component  systems  described  above.  However,  the  need  to  control  two  adsorbate 
coverages  and  the  large  range  of  relative  coverage  possibilities  makes  such  studies 
difficult. 

Lad  et  al.  [84]  have  studied  0  and  5  adsorbed  on  Ni(  1 1 1 )  to  a  total  coverage  of 
one-half.  A  mixed-phase  coexistence  region  with  O-rich  and  5-rich  c(2x2)  domains 
forms  at  low  temperatures  and  transforms  into  a  c(2x2)  phase  with  substitutional 
disorder  as  T  is  increased.  At  higher  temperatures,  the  system  is  disordered. 

5.  SYSTEMATIC  BEHAVIOR 


Enough  phase  diagrams  have  been  determined  to  begin  examining  the  systematic 
variations  between  them.  This  offers  the  possibility  of  much  more  stringent  tests  of 
theoretical  models  than  is  possible  in  an  in-depth  study  of  a  single  adsorbate- 
substrate  system. 

The  chalcogens  0,  S,  Se  and  Te  adsorbed  on  Ni(  111)  have  been  studied  in  some 
detail  [69,84,88],  The  substrate  bonding  is  believed  to  be  similar  in  each  case  but  the 
adsorbate-adsorbate  interactions  will  vary.  Figure  7  shows  the  partial  phase 
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diagrams  for  0,  5  and  Te.  Surprisingly,  these  diagrams  differ  significantly.  For 
example,  a  p(2x2)  phase  is  the  ordered  phase  at  lowest  coverage  for  0  and  5  whereas 
Te  forms  a(2/3x  2/3)  phase.  A(/3  x  /3)  phase  occurs  below  room  temperature  for 
0,  is  not  observed  for  S  and  exists  only  above  700  K  for  Te.  Complex  coincidence 
structures  form  at  high  coverages  for  5  and  Te  while  0  forms  bulk  compounds.  The 
precise  origins  of  these  complex  differences  is  not  understood  and  offer  new 
challenges  to  theory. 

6.  CONCLUSIONS. 

This  review  has  attempted  to  summarize  recent  experimental  studies  of  critical 
behavior  in  continuous  order-disorder  phase  transitions  on  surfaces.  Several  points 
should  be  emphasized  in  conclusion. 

First,  the  appropriate  proceedures  for  acquisition  and  analysis  of  the  diffraction 
data  are  only  beginning  to  be  understood.  Surface  preparation  methods  are  relatively 
crude  and  important  contributions  of  finite  size  effects  introduced  during  the  sample 
orientation,  polishing  and  cleaning  cannot  be  ruled  out  in  any  of  the  systems  studied  so 
far.  As  demonstrated  in  the  cases  of  AuC  1 10)  and  W(001)  the  effects  on  the  critical 
behavior  can  be  significant.  This  may  be  the  most  important  problem  facing 
experimentalists  interested  in  critical  phenomena  on  surfaces. 

More  accurate  methods  are  needed  for  determination  of  absolute  coverages.  None 
of  the  methods  employed  have  absolute  accuracies  better  than  10-20  percent. 

LEED  has  been  the  only  experimental  technique  used  for  quantitative  studies  of 
critical  exponents.  In  addition  to  the  problems  associated  with  correction  for 
background  and  thermal  diffuse  scattering,  which  will  be  common  for  any  diffraction 
method,  the  important  question  of  whether  multiple  scattering  changes  significantly 
near  a  phase  transition  remains  unanswered.  Combined  LEED  and  Glancing  Incidence 
X-ray  Diffraction  measurements  carried  out  on  the  same  system  should  clarify  this 
point  in  the  near  future. 

The  method  of  Bartelt  et  al.  [15]  which  uses  averaged  LEED  intensities  to 
determine  the  exponent  a  is  not  sensitive  to  multiple  scattering.  Furthermore,  the 
necessary  measurements  can  be  carried  out  with  a  low  resolution  LEED  apparatus. 

This  is  an  important  development.  It  means  that  even  the  poorest  quality 
instrumentation  can  be  used  to  study  some  aspects  of  critical  behavior  in  a 
quantitative  way. 

Finite  size  effects  and/or  instrument  quality  have  so  far  limited  the  nearness  to 
which  the  critical  temperature  can  be  approached  to  a  few  parts  per  thousand.  This 
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means  that  higher  order  terms  in  the  expansions  of  Equations  (8H12),  which  define 
the  critical  exponents,  might  not  be  ignorable. 

Finally,  it  is  curious  that  nearly  all  of  the  surface  systems  studied  so  far,  Table 
2,  have  critical  exponents  which  are  most  consistent  with  those  of  the  Ising 
universality  class.  Clearly,  the  hope  of  finding  real  physical  examples  of  all 
universality  classes  has  yet  to  be  realized.  Although  progress  toward  using  surface 
phase  transitions  as  models  of  two-dimensional  critical  phenomena  has  been 
substantial,  numerous  problems  remain  and  much  work  needs  to  be  done  before  this 
important  subject  is  placed  on  a  firm  footing  solidly  based  on  experimental  evidence. 

In  spite  of  these  current  difficulties,  recent  progress  has  been  substantia!. 
Improved  techniques  will  develop  rapidly  and  experimental  studies  of  phase 
transitions  are  certain  to  have  an  exciting  future. 
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Figure  Captions: 


Figure  1:  Thermodynamic  phase  diagram  for  submonolayers  of  selenium 

adsorbed  on  a  Ni(00 1 )  surface  [5].  The  open  dots  are  experimentally 
determined  phase  boundaries.  The  solid  lines  and  vertical  dashes 
mark  the  phase  boundaries.  All  phase  boundaries  are  believed  to  be 
continuous  for  this  system. 


Figure  2:  Structural  models  of  the  ordered  phases  of  Figure  1.  The  solid  dots 
represent  Ni  atoms  in  the  surface  layer  of  the  (001 )  substrate.  The 
circles  represent  Se  atoms  adsorbed  at  four-fold  symmetric 
hollows  on  the  substrate  lattice. 

a)  The  p(2x2)  lattice. 

b)  The  c(2x2)  lattice. 

c)  Selenium-selenium  interaction  energies.  The  phase  diagram  in 
Figure  1  requires  E3  <  E2  <  Ej. 


Figure  3:  a)  A  (2x1)  structure  on  a  square  lattice  substrate.  Two  symmetry 
equivalent  domains  are  shown;  their  boundary  is  marked  by  the 
dashed  line. 

b)  Surface  Brillouin  zone  for  (a).  The  substrate  reciprocal  mesh 
vectors  are  aj*  and  a2*.  The  zone  center  is  at  r.  The  reciprocal 
vectors  k|  and  k2  describe  the  order  parameters  for  the  two 
(2x1)  domains.  The  wavevector  q  describes  the  order  parameter 
of  a  c(2x2)  structure  as  shown  in  Figure  2(b). 


Figure  4:  Temperature  dependence  of  the  components  of  the  diffracted 

intensity  near  a  continuous  order-disorder  phase  transition.  The 
upper  portion  shows  the  intensity  due  to  long-range-order  I|_rq  and 
the  amplitude  of  the  susceptibility  x(OJ),  both  measured  at  the 
Brillouin  zone  center.  The  lower  portion  shows  the  angular  profiles 
for  temperatures  just  below,  at  and  just  above  Tc.  Only  single 
scattering  is  assumed. 
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Figure  5:  Critical  behavior  of  the  Au(l  10X1x2)  order-disorder  phase 

transition  [  1 7],  (a)  The  temperature  variation  of  the  components  of 
the  intensity  due  to  long-range-order  (•)  and  to  fluctuations  (+). 
The  solid  line  is  the  behavior  predicted  by  Onsager  for  a  simple 
Ising  model  [53].  (b)  Log-log  plots  of  (LR0  (•),  x^,T)  (+)  and 
ASj|(^)  (x)  vs  reduced  temperature.  The  slopes  yield  the  critical 
exponents  13,  y  and  v  respectively. 

Figure  6:  Schematic  phase  diagram  for  the  XY  model  with  cubic  anisotropy. 
The  parameter  w  measures  the  strength  of  the  anisotropy,  w 
positive  corresponds  to  a  (2x2)  or  (2x1 )  ordered  phase  and  w 
negative  to  a  c(2x2)  ordered  phase.  Along  the  order-disorder 
transition  line,  J3  varies  continuously.  (After  reference  61 ). 

Figure  7:  l  xperimental  phase  diagrams  for  0,  S  and  Te  adsorbed  on  Ni(  1 11 ). 
Dashed  lines  are  first-order  transitions  and  solid  lines  are 
continuous  transitions.  (After  references  69,  84  and  88). 


TABLE  2 


SYSTEM  a 

J3 

y 

V 

TC(K) 

References 

Audi  0)(  1x2)  0+0.2 

0. 1 3i0.02 

1.75+0.03 

1.02+0.02 

650-713 

16,35 

Si(1 1  0(7x7) 

0.22+0.02 

* 

* 

1  140 

47 

0-Ni(l  1 1)p(2x2) 

0.14+0.02 

1.9  +  0.2 

0.94+0.10 

All 

68,69 

0-W(112)p(2x2) 

0.13+0.01 

1.79+0.14 

1.09+0.1 1 

900 

80 

H-W(l  10)(2x1) 

0.13+0.06 

78 

H-W(1 1 0X2x2) 

0.25+0.07 

78 

Ba-W(l  10X3x2)  «0+ 

125 

81 

Cl-Agd00)c(2x2) 

0. 1 2+0.04* ** 

77 

*No  critical  scattering  was  observed. 

**Quantity  actually  determined  is  13/d-ct). 
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